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Abstract 

Recently it was suggested that the problem of species doubling with Kogut-Susskind 
lattice fermions entails, at finite chemical potential, a confusion of particles with an- 
tiparticles. What happens instead is that the familiar correspondence of positive-energy 
spinors to particles, and of negative-energy spinors to antiparticles, ceases to hold for the 
Kogut-Susskind time derivative. To show this we highlight the role of the spinorial "en- 
ergy" in the Osterwalder-Schrader reconstruction of the Fock space of non-interacting 
lattice fermions at zero temperature and nonzero chemical potential. We consider Kogut- 
Susskind fermions and, for comparison, fermions with an asymmetric one-step time 
derivative. 



1 Introduction 



The implementation of finite baryon density on the lattice [|], ||] has proven to be a 
hard task. Several proposals have been made to explain the unexpectedly low chemical 
potential fx at which the baryon density and other quantities start to grow. The onset 
at a // much lower than the value ton/3, as expected for a free nucleon gas, could be due 
to the enhanced number of flavours (when using Kogut-Susskind fermions) producing a 
much stronger nuclear binding than in nature 0, || . 

As the onset is so early that its energy can be almost confounded with the existence of 
a 'baryonic' Goldstone mode, it was also interesting to search for possible enhancement 
of mesonic modes due to some kind of lattice artifact. In |J it was pointed out that 
by selecting only the positive-energy poles of the free Kogut-Susskind propagator one 
still obtains energies and baryon number densities composed of particles as well as 
antiparticles. This is in contrast to Dirac fermions in space-time continuum where 
positive-energy spinors always correspond to particles and negative-energy spinors to 
antiparticles. As observed in ||, that feature of continuous fermions is reproduced on 
the lattice if an asymmetric, one-step time derivative is used. (In fact, such fermions 
have the same Hilbert space as Susskind's fermions || on the Hamiltonian lattice.) 
It was then suggested in || that the Kogut-Susskind action does not allow a clear 
distinction of particles and antiparticles. That conclusion, however, is not justified. In 
fact, particles and antiparticles can be reliably identified only from the Hamiltonian and 
baryon number operators in Fock space which we therefore reconstruct below. Apart 
from the usual doubling of flavours, nothing peculiar will be seen here to occur with 
Kogut-Susskind fermions. Rather, it will be the correspondence of positive/negative- 
energy spinors to particles/antiparticles which is destroyed by the Kogut-Susskind two- 
step time derivative. 

Intuitively, in this respect, the essence of the technical, field-theoretic arguments 
presented in the next sections is as follows. In Euclidean time r, "particle" solutions 
of the Dirac equation are those which remain bounded for r — > +00 while "antipar- 
ticles" remain bounded for r — > —00; this corresponds to the Minkowskian notion of 
particles moving forward and antiparticles backward in time. The rule is implicit in 
the Osterwalder-Schrader reconstruction algorithm (cf. remark following eqn. (|l6|)). For 
the Dirac equation in continuous Euclidean time, ip = —Hip, this implies that the pos- 
itive eigenvalues of H are associated with particles and the negative eigenvalues with 
antiparticles. Now, if time is discretized, the Dirac equation is converted into a recur- 
sion relation. The most straightforward, asymmetric discretization would give (with ax 
being the lattice spacing and the spinorial Hamiltonian diagonalized) 

Vv+i -ipr = -a T E tp T 

Depending on the size of OT-Ecutoflf this will distort the continuous-time dependencies 
somewhat but it will preserve the qualitative features. In contrast, the symmetric dis- 
cretization would give the recursion relation 

i' T+ i - tJj t -! = -2a T E tp T 

which is for any ot a qualitatively different type of equation. In fact, looking for 
"particle" solutions with an exponential decay e _tJr in the positive time direction, we 
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now find two of them: 




,..,,_.= + sinhu; > 
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A similar degeneracy, involving both signs of spinorial energy, is obtained with solu- 
tions bounded for r — ► — oo. Hence, positive spinorial energies are no longer reserved 
for particles, but occur as antiparticle solutions as well. Of course, those arguments 
require formal substantiation. In the field theory that we consider here, particles and 
antiparticles will be clearly identified by the baryon number. 

The Hilbert spaces corresponding to space-time lattice fermions have all been recon- 
structed long ago as part of the issue of the positivity of the transfer matrices. This 
was done in |?], [8L || for Wilson fermions, and in [ 10 1 for staggered fermions. In all 



those references, a direct comparison was made between matrix elements in a standard 
fermion Fock space and the Grassmann functional integral. While this method enables 
one to deal with interacting fermions from the outset, it does not assign a Fock-space 
interpretation to the path-integral variables themselves. For Kogut-Susskind fermions, 
in particular, auxiliary Grassmann integrations had to be introduced, preventing the 
Fock space interpretation of the original variables. In the present context, where we 
want to identify the role of positive- and negative-energy spinors (the Fourier trans- 
forms of the path-integral variables) in the fermion Fock space, we therefore prefer to 
work our way through the Osterwalder-Schrader reconstruction scheme [[ll], |l2| which 
does not anticipate any knowledge of the Hilbert space. For a detailed introduction to 
this formalism we refer the reader to we shall only need the simplest of the formal 
concepts, and the mathematical tools required in the following will not go beyond some 
suitably organized algebra of Grassmann derivatives. 

Fortunately, the present issue allows us to restrict ourselves to non-interacting fermions 
at zero temperature, but at finite chemical potential /i. 

Furthermore, we will restrict our examples to cases where the inverse fermion matrix 
can be given explicitly, in terms of spinorial energy eigenstates and eigenvalues. Thus 
we consider fermion actions of the general form 

S = tp n K nn ftp nl n = (n 1 ,n 2 ,n 3 ,n4) = (n,n 4 ) (1) 

where K nn > = Sn,n'D((j,) n n , + S n „/ — E HH , 

4 ' 4 os 

such that the time-changing (derivative) part D commutes with the equal-time, spinorial 
"energy" part E. This form can be achieved both for Kogut-Susskind and for time- 
asymmetric || fermions. For Wilson fermions, such a simplification does not seem to 
be possible. 

For the asymmetric fermions it will be advantageous to allow for an anisotropy of 
the space-time lattice; hence the appearance in ([!]) of the spatial and temporal lattice 
spacings, as and ax- 

In section 2 we recall the basic steps in the procedure of Osterwalder-Schrader re- 
construction, and the standard tools we have at hand for non-interacting fermions. In 
section 3 we apply this, in some detail, to fermions with a one-step time derivative. 
We thus recover the familiar zero-temperature scenario at finite chemical potential. In 
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section 4 we provide the analogous expressions for Kogut-Susskind fermions. A com- 
parison of the operator expressions for the energy and for the baryon number in Fock 
space (as opposed to the linear space spanned by the Grassmann variables) will show 
that a chemical potential enhancing the baryons will deplete the antibaryons. In the 
concluding section 5 we discuss relations with the positivity of the transfer matrices and 
with finite temperature. 



2 Basic ideas of reconstructing H and B 

From the path-integal expectation values and correlations we wish to recover the multi- 
fermion Hilbert space, the scalar product of the Hilbert space vectors, and the operators 



of time evolution and baryon number. The basic ideas are as follows [13]. Both bra 
and ket vectors contribute to (• • -) P ath integral- The ket vectors are the functionals of 
Grassmann variables ip T and ip T with Euclidean times r > (more precisely, they are 
equivalence classes of such functionals; cf. below). The bra vectors are the same with 
r < 0. Thus bra and ket functionals are related through a time-reversal operation 0. 
On Grassmann variables it acts as 

e^(r) = $(-t) 0$(r) = Tp(-r) 

and is extended as an antilinear operation on functionals of ip and tp. We shall adopt 
the rule that complex conjugation of a product of Grassmann variables implies reversed 
ordering of factors. 

The scalar product of a bra and ket vector is given by a path-integral correlation: 
if F and G are ket functionals, and \F) and \G) their interpretations as Hilbert space 
vectors, their scalar product is 

(F\G) = (@F G% ath intcgral = ^J [d#d$ e s 0F[^V] T <o G^VUo 

For a quantum-statistical interpretation the scalar product must satisfy (F\F) > 
which it does in all the cases we consider. Null-norm functionals (which have (F\F) = 
although F ^ 0) can then be factored away. This will, in fact, simplify the structure of 
the Hilbert space. 

We will extract the transfer matrix operator T from the path integral correlation 
(F\T\G) in such a way that it will reproduce these matrix elements for arbitrary F and 
G. The transfer matrix, which defines the true energy (as opposed to spinorial "energy" ) 
by T = e~ H , derives from a unit shift of the time index on all the Grassmann variables 
in a ket functional 

( FfaM I f | G[Vv,?/v] ) = ( F[</v,W] I G[Vv+i,W+i] ) 

The baryon number operator will likewise be identified from its matrix elements (F\B\G) , 
i.e. from the path integral correlation (OF B G) in which B[ip,ip] = dS/dfi is the usual 
observable. 
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It will be essential in the following to have an explicit representation of Grass- 
mann-integral expectation values and correlations at hand. In the most general finite- 
dimensional case, if ip n and xf) n are the integration variables and if 0[if),t/j] is any ob- 



servable or correlation in terms of them, then [14] 



integral a \- 



inn 



exp 



Yk- 1 



d d 



\n,m 



(2) 



' dijj n dt/j n 

To keep the notation as simple as possible, we shall refer only to Grassmann derivatives 
"from left to right" {d/dip = <9). The reader who wishes to check the crucial signs for 
himself will find it convenient to use Grassmann derivatives "from right to left" as well, 
noting that d d FG = -Fdd G 

We shall be only concerned here with the Euclidean time structure of the inverse 
action kernel K~^, where n = (ni, ri2, 77,3, 77.4) etc. We shall use two different labelings 
of time slices: n& = 0, ±1, ±2, ... as appropriate for the space-time functional-integral 
approach, and r = ±1, ±2, . . . which will emphasize later on the correspondence between 
bra and ket functionals. For later convenience we note that 



T = 714 for 77-4 > T = 724 — 1 for 714 < 

With respect to the spatial indices, we transform into a basis in which the equal-time 
part of the action is diagonal. For Kogut-Susskind fermions this involves spatial Fourier 
transformation and linear combination of the momenta at the edges of the Brillouin 
zone. Let p denote a suitable (multi)index for the basis vectors. For example, p can be 
chosen as 



P = (Pi,P2,P3,cr,s,f) Pk& 



7T 



7T 



2a s ' 2a s 



±1, 



4 f 



u,d (3) 



where p}~ are the components of spatial momentum. The rest of the Brillouin zone, 
reachable by addition of "edge" momenta (±7r/a, ±7r/o, ±7r/a) corresponds to the de- 
grees of freedom of spinorial "energy" sign a, spin s, and flavour / |lE}| . In the following, 
<j(p) will be of particular interest to us. Using 



E(p) = sin 2 p fe + m 2 



the spinorial "energy" is a(p)E(p). 

We can now rewrite the action kernel of equation @) as 

Knn> = 8 P , P > yD(n) n4 ^ 4 +5 n4in > 4 —a(p)E(p)j 

In the following we reconstruct the fermion Fock space, and the operators of energy and 
baryon number in particular, for two versions of the time derivative. 
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3 Asymmetric time derivative 



We first analyse fermions with a one-step, asymmetric time derivative. This will lead to 
expectation values (and operators) for the energy and baryon number resembling those 
familiar from space-time continuum [||. 

The time-asymmetric action kernel, corresponding to equation (2.15) of ref. [n], is 

la 3 

K nn , = ^mS nn > + - f^ £ 7fc (S r +k - 6 ,_ h ) + 74 {^b n n , +l - 5 n y) 
as * a s k=1 

where the Euclidean gamma matrices are characterized by {7^, 7^'} = 5 VU > and 7J = j v . 
To reduce the fermion degeneracy, the 4x4 matrix structure can be eliminated by first 
absorbing a 74 factor into the tp integration variables and then diagonalizing in spin and 
restricting to one of the components by analogy to Kawamoto-Smit. Hence let us define 
«fc = ililk and 

X(n) = ${n) 74 ar«2 2 «r x(n) = a^a^af^n) 
With respect to the x an d X spinors the action kernel is 

la 3 

K nn > = ^m5 nn /7 4 r 4 (n) + £ T k (n) (6 , +i - 5 ,_ k ) + ePS l+i - 5 n y 
as * a s k=1 

where T v (n) = (— l)™ 1 " 1 ^-i. A 74 matrix has remained in the mass term but, as- 
suming it to be diagonalized, the spinor components are now decoupled so that we can 
build a flavour-reduced fermionic system by projecting onto the first component. We 
note that the equal-time part of this action coincides with Susskind's spatial fermions 
||. After spatial Fourier transformation the asymmetric kernel becomes 

K nn > = S pp/ (e**tf n4+1)Tl £ - 5 n4: n'J + a(jp)E(jp) — 5 n4 y 4 I 

To evaluate any scalar product or matrix element, such as ((@F)BG), using (^), we 
need to know the inverse of the action kernel. Using 

e = In (l - ^L a (v)E(v) 
V as 

we obtain 

\ -e-^e^-^^-O ifn 4 <n 4 K , ... 
K ™' = \ ifnI><) V iOTe> " ^ 

x _ \ if n 4 < n' A K , , 

K nn' - I e -e e (6- M )(n 4 -<) if „ 4 > ^ j V tor e < l 5 J 

It should be mentioned here that the zeros above will produce zero-norm states later 
on. Right at e = fi the inverse kernel is unambiguously determined by the antiper iodic 
boundary condition with respect to time. The appropriate limit is 
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It is convenient to adapt operation (Q) to the special case of O = F-G with F a negative- 
time and G a positive-time functional. (This also covers the case of a matrix element 
like (@F B G) .) Sorting out the derivatives acting on F or G exclusively, we can write 

»pE^^ F =M« +G )) (6) 
where the C stand for "contractions" and 



v t,t' >0 
t, t' < 



<9 W C^pr' 



while the remaining "mixing" term is 



Co = exp £ (-e~n ■ 

p: <:>/-* \t>0 



p: £>M \r>0 °W P -tJ \ , >q GTp.r' 




The operations C+ and C_ only modify the ket and bra functionals, respectively. The 
operation intertwining kets and bras is Co, and the last step is to discard all remaining 
■0s and ?/>s. It will be most convenient to formulate the fermion Hilbert space in terms 
of F' = C-F and G' = C + G, i.e., at the Co stage after the C+ and C_ operations. 

We are aiming at operator expressions for the energy (transfer matrix T or Hamilto- 
nian H) and for the baryon number. As C+ and C_ depend onr-r' only, they commute 
with the time evolution. Thus T can be determined at the Co stage immediately; see 
below. The baryon number operator is most easily identified from the path-integral 
observable "at r = 0" (i.e., involving r = — 1 and r = 1) which is 

OS - 

B = j- = £ e^n,-i^n,i = (8) 

°^ n V 

In the second equality we changed to the eigenbasis of the spinorial energy, and arranged 
Grassmann variables in time order. V> P) -i and ip Pt i can now be interpreted as multiplica- 
tion operators acting on bra and ket functionals, respectively. Thus the matrix elements 
of the baryon number operator are 

(F\B\G) = e" <0(^,i.F) V P ,iG> P ath integral = e"C {C_ {Q^ Ptl F)C^ Ptl G} 

To bring the field operators (multiplication operators) of the above equation to the 
Co stage, we apply the Baker-Hausdorff formula, obtaining for example the following 
primed form (with any ket functional G) 

w =L l + {l : >;;L e -,v><-^W (9) 
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On the bra side, we obtain the G-reflected expression. In a similar way, by commuting 
ipa and v/js through Co one obtains the adjoints of field operators. For example, we have 

d \ 




r>0 Wp>tJ 

The ijjp—i on the RHS is discarded when all ij) and V> are set to zero in the final step 
of the scalar product. The bracketed expression on the RHS is the ket-adjoint to the 
bra-multiplication by ipp —i- To further develop our example, this can be used to rewrite 
Co {{Qip p ^F')ijj p ^G'} as Co ^{QF')il^ p l ip p ^G'^ which is part of the baryon number oper- 
ator we are interested in. From an expression such as the last one, in which all operator 
action is thrown on the kets, one finally abstracts away the functionals and the scalar 
product operation. 

On the Co stage, as we see from (0) and @, Grassmann derivatives occur only in 
those linear combinations of time indices which, for each momentum p, correspond to 
"physical" time dependency. Linear combinations of ip T and ip T orthogonal to it drop 
out from all matrix elements (they create the null-norm functionals). The only relevant 
variables are oc J2t>o e^~^ T and *ff(p) oc J2t>o e ^~ e ^ T ■ They are conveniently 

defined by 

V e (M-E)r g = , g 



T>0 



so that, in the ket functionals, we have 

4>p,t = e^ _e > T tf (p) + orth. e> fi . . 

if, PtT = e ^-^ T ^(p) + orth. e</i [ ' 

(The and ^(p) actually represent equivalence classes.) Omitting the irrelevant 

"orthogonal" parts we obtain functionals of ^(p) and *f>(p) only. We note that the ^ 
are non-null only for e < /U (in ket functionals) and the ^ only for e > /x. This will 
be different for Kogut-Susskind fermions, leading to an additional doubling of species 
there. 

In terms of the new variables, the baryon number (||) for time-asymmetric fermions 
takes the simple operator form 

At the Co stage, the transfer matrix is also easily identified. If all Grassmann variables 
in a ket functional are shifted by one time step, equation ([[(]) shows that ^f(p) picks up 
a factor of e^~ e , and ^f(p) a factor of e e ~^. The generator of such a transformation is 

h — fj,B = — log r = £ ( £ - mp)^ + £ (-e + m)*(p)^) (12) 
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We can rewrite this in the standard Fock-space form, using e m = log(l + ^-m): 

H-^B = (e-n)b*(p)b(p)+ £ (e + f ,)b*(p)b(p)+ £ (e - M ) (13) 

p:e>e m p:e>t m p: p>e>e m 

B= £ (6*(p)6(p)-6*(p)6(p))+ E 1 ( 14 ) 

p:e>e m p:e<-e m 

where we have used e m = log(l + ^m) and defined creation and annihilation operators 
for antibaryons (denoting by p the multiindex p as in (||) with a reversed sign of spinorial 
energy) by 

6*(p) = *(p) 6(p) = 5/a*(p) e(p)>e m 

and for baryons by 

W \ 0/0* (p) a* > e > e m ° W \ *(p) /x > e > 6 m 

The 6* and b* satisfy canonical anticommutation relations with b and b, implying stan- 
dard identification for number operators. However, for e(p) or fi large (comparable to 
the uv cutoff) they deviate from the hermitian adjoints of b and b by a real factor 

P(p) = e - 2t -%*{p) foralle>e m 

b\p) = e' M b*(p) e>n tf(p) = e> M - 2e b*(p) /i > e > e m 



As we see from ( |l2| ) the functional representing the ground state at chemical potential 
fjL is Go = 1- in terms of Fock space operators it is characterized by 

b(p)\fi)=Q allp 

b(p)\n)=0 fore>/i (15) 
tf(p)\lA = for fj, > e > e m 

Thus the ground state is filled with baryons of all momenta, spins and flavours with 
energies in the range fx > e(p) > e m . 

The zero-point energy in ( |l~3| ) is due to the above definition of H — [iB as a plain gen- 
erator of a temporal translation. Had we defined the baryon number operator as a plain 
generator of phase rotations, we would have avoided the divergent but ^-independent 



contribution in (14) which we must now subtract by hand. 

It should be noted that definition (|i~3|), based on a single time step, is valid only 
if e in the above definition e = ln(l — ^-a(p)E(p)) is a real number for all energy- 
momentum indices p. Since the maximal eigenvalue of the free staggered Hamiltonian is 
V3 + m 2 , m being the fermion mass, the condition is ax < as/\/3 + m 2 . If the temporal 
spacing is too coarse to satisfy this, then one can still define a two-step Hamiltonian, 
H — fxB = — g logT 2 as customary for staggered fermions. 
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4 Kogut-Susskind time derivative 



We now consider the two-step time derivative as it is used with Kogut-Susskind fermions 
|[], Q]. The action kernel, allowing for a space-time anisotropy, is 

1 3 

K nn , = ^m5 nn , + T k( n )( S n +k,n' ~ S n -k,n') + T ^ n )^ 5 n+ln> ~ e ~^n-4,n') 

«S «S fc=1 

Absorbing a factor of r 4 into the ip integration and Fourier transforming with respect 
to n we arrive at 



K nn > = S pp/ K ni n' 4 {p) where K nin ' 4 (p) = I {e^5 nA+ iy A - e ^-i,^) +a(p)E(p)^5, 
Let us restrict to /j, > and put 



a s w n 4 n' 4 



arsinh ( — oE 



The inverse kernel K 1 , (p) then reads 

714,714 « ' 



-M(n4-n 4 ) e - e |i4-n 4 l f 1 for n 4 > W 4 

_ X 1 (-l)^-" 4 for n 4 < n' 4 6 > M 



coshe 



e -u( ni -n' 4 ) e e\n 4 -n' 4 \ ( ^^r^-^+l for ^ > R , 



coshe ) —1 for n 4 < n 4 



e < -// (16) 



-u( n4 -n>)e~~ <n4 ~ n ' d - {-l)^-<e< n ^ f 1 for n 4 > n 4 
cosh e I for n 4 < n 4 

A general remark may be in order here. The inverse kernel consists of positive-time and 
negative-time solutions of the discretized Dirac equation, suitably patched together at 
n 4 = n' 4 . In eqn. (|2|) the index n 4 always goes with a d/dip. If for a particular solution 
re 4 runs over the positive half-line, then ^s are detected in the ket functional; if n 4 runs 
over the negative half-line, the ips are in the bra functional where they are interpreted as 
complex-conjugates. Whatever corresponds to a "particle" in one case thus corresponds 
to an "antiparticle" in the other. 

To carry on with the technicalities, we change to time labels r, t' and split the sum- 
mation over time according to (||). In particular, the intertwining operation analogous 
to (0) now reads 



cosh e \ dib v T I I , n dib„ _ r > 

p: £>-/! \r>0 / \r'>0 
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From this we can identify the relevant Grassmann variables at the Co stage of the scalar 
product. The procedure is largely the same as with asymmetric fermions. Again, the 
new variables can be conveniently defined by 

\ e+fl \ T _d_ 8 \- -|e+ M lr 9 _. 9 



Y e ~ |£+/j|r =• - V e 



so that, in the ket functionals, the Grassmann variables decompose as follows: 

e > n i> PiT = (-l) r e^ -e ) T *(p, -) + orth. ^ PiT = e~ (e+M)r #(p, +) + orth. 
|e| < /U ^ P ,r = + orth. V P>T = e-<- e+ ri T V(p, +) + (-l) T e< e -") T ¥(p, -) 

+ orth. (17) 
e < -ii V5 PiT = e( e+ ^ T *(p, +) + orth. ip P)T = (-1) T 'e( 6 ~^ T ¥(p, -) + orth. 

It is interesting to note that the limit ot — ► can be taken only after these trans- 
formations when the (— 1) T time-dependency mode has "transmuted" into an internal 
degree of freedom. The same kind of combinations have to be considered in the spatial 
directions giving their flavour degeneracy. However, as one is ultimately interested in 
interacting fermions which allow the continuum limit to be taken only at the very end 
of the calculations, the alternating modes of free fermions do not a priori signal any 
difficulties. 

The baryon number as a path-integral observable |l|, [2| on the zero-time slice is 

dS 

B = TP = ( e/ V P ,-iW,i + e-'Vp.i^p.-i) 
a P p 

An intermediate result for the baryon number operator at Cq level is 

- (-db<^ + *<* + >) (- e ^ (p '- )+ ^^ry) 

- £ (*(p.+)+*(p.-))f^«^T+ r "" " ' 



2 77' \ coshe d^f(p, +) cosh e d^f(p, — ) I 

p: \e\<fi \ f 

+ 2 eAt p: ^("coshe^(p,-) +*^' + )) ("^ M * (p ' _) + cosheO*(p,+) 
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I -m V" ( 1 9 q,( _\\ ( ^+^( e^~ £ d \ 

2 6 p: ^Uoshe5^(p,+) + (P ' 7 V ^' +) + coshed*(p,-)J 

This actually simplifies to 

+ ^O-^-'^)-^^'^)) (18) 

The transfer matrix of Kogut-Susskind fermions is known to be not of the preferable 

. Nevertheless one can define a quantum field theory 
using the two-step Hamiltonian H = — \ log T 2 . This is hermitian whenever T is. Thus 
we obtain, generalizing to /i > 0, 

* - " 6 = P £ ( (£ + " ,<p(p ' + (£ - -><»&=)) 

+ ( ( " £ + _) ^) + ( ~ £ " fc+) aS^+o) 

This can be rewritten in Fock form by appropriate identification of baryon and an- 
tibaryon operators. Let us put e m = arsinh (^-m) and rj = db. Furthermore, to indicate 
the relation between timelike "flavour" and the sign of the spinorial energy, let rjp and 
— r\p stand for (pi,P2,P3, ±c, s, /) and (pi,P2,P3, T ") s, /), respectively. The antibaryon 
annihilation and creation operators are 

b(p,r)) = d/d^{-qp,r}) b*(p,rj) = 9{r]p,r]) e(p) > e m 

The corresponding baryon operators are 

b{p,rj) = d/d^f(—rjp,ri) b*(p,7]) = ^f(—rjp,rj) e(p) > \i 

b(p,rj) = W(-r)p,ri) b *(P, v) = d/dW(-T]p, rj) p, > e{p) > e m 
We then have 

B = Yl X) ( b *(P>V) b (PiV) ~ b*(p,rf)b(p,r))) 

p: e>e m T) 

and 

H-/jB = T,(( e -^(P^)b(p^) + (e + p)b*(p,7 ] )b(p,7 ] ))+ Yl 2 ( e -^) 

p:t>t m V p:^>e>e m 

We note that creation and annihilation operators are only defined for multi-indices p 
with a positive energy. In terms of these operators the Hamiltonian and baryon number 
take the standard form of baryons and antibaryons with an extra two-component flavour 



form e with H hermitian [10 
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rj. The chemical potential ji > decreases the energy of baryons and increases it for 
antibaryons, as expected. 

As in the asymmetric case, the b* and b* at large e or /i (comparable to the uv 
cutoff) deviate from the hermitian adjoints of b and b: 

Hp, rj) = e _e_A1 cosh e b*(p, rj) e > e m 

I r(v,rj) = e~ e+Ai coshe b*(p,rj) e > fi b'(p,r]) = — b*(p,rj) fi > e > e m 

cosh e 

where e = e(p). The hermiticity of H is not affected by those relations. 

As with asymmetric fermions, the ground state functional is Go[\l/,^] = 1. In 



terms of creation and annihilation operators this implies the same relations as in (15) 
for both of the time-like flavours. The ground state has a nonzero baryon density if 
sinh// > ^rn. Thus, again, fi and m are related in a nonlinear way, reminding us of a 



degree of ambiguity in the quantification of chemical potentials on the lattice 16, 17] 



5 Conclusions 

We have elaborated on how the sign of the "energy" of Grassmann functional integral 
variables is associated with the Fock-space degrees of freedom in the two cases of Kogut- 
Susskind and time-asymmetric lattice fermions; the central equations in this respect are 
( [l~0[ ) and (^). Our main conclusion is that the usual implementation of a chemical 
potential for the baryon number, namely through real-valued propagation factors in the 
discretized time derivative for any kind of fermions, will indeed enhance or deplete the 
baryons and antibaryons just as it is supposed to do. 

It is quite plausible that time evolution with Kogut-Susskind fermions should involve 
some peculiarities. In continuous Euclidean time, the positivity of the transfer matrix 
is guaranteed already by the basic requirement (F\F) > [|ll], §. By contrast, the 
Kogut-Susskind one-step transfer matrix has both positive and negative eigenvalues at 
any time-like lattice spacing [FLCl] and is therefore usually regarded as involving both 



time-evolution and flavour rotation [15]. In fact, the hermiticity of the energy in this 
case is accomplished by simply discarding the sign of the negative eigenvalues of T, 
namely in defining H = — ^logT 2 instead of H = — logT. The fact that this leads 



to a sound quantum field theory in principle has been well-known [10|. In the practice 
of Monte Carlo simulations, when dealing with propagators of particles composed of 
Kogut-Susskind fermions, a contribution of "odd parity" states (involving (— 1) T ) is 
therefore fitted to a free Kogut-Susskind particle with the same flavour rotation and 
time evolution. 

Complications arising from the nonpositivity of the transfer matrix are avoided 
with Wilson fermions J|], and also with the asymmetric fermions of section 3, pro- 
vided that the temporal lattice spacing is sufficiently smaller than the spatial one. For 
non-interacting fermions the condition is ax < as/V3 + m 2 ; similar bounds exist for 
fermions interacting with a compact gauge field. 

In computational practice, and related analytical study such as [||, |l(], ||], the en- 
ergy expectation value e of lattice fermions at finite temperature is not derived from 
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the logarithm of the transfer matrix, e = tr He~~ HN/3 , but by differentiating the ac- 
tion with respect to the temporal lattice spacing. In terms of Fock space operators 
this can be obtained from e = tr ET N > 3 = tr Ee~ HN ' 3 where E is the operator recon- 
structed from the path-integral observable J2 P ^,1-^(^)^,1 f° r one-step fermions, or 
\ J2 P {$p,-\E{p)nl) p -i + , i/jp^E(p)ipp j i) for Kogut-Susskind fermions. In fact, E and H 
are different only for momenta comparable to the uv cutoff so that the use of E to 
obtain an energy density is certainly justified by its computational advantages. 

It is reassuring to find that particles and antiparticles are safely recognized by their 
internal transformation properties (of which we have considered the baryon number) and 
not by any kind of spinorial energy. The latter would be an ill-defined concept, indeed, 
in the presence of dynamical gauge fields. The Hamiltonian of interacting fermions can 
only be reconstructed as part of the complete, translation-invariant fermion-boson field 
theory. Naturally, such a reconstruction had to emphasize local terms (as it did in 
[0, ||, 10 1) rather than global diagonalizations. 



As the time-asymmetric action could provide an exact algorithm with two flavours, 
it will be interesting to study it further in the SU(3) interactive case, in concordance 
with our view that the early onset problem could be due to the excessive proliferation 
of degenerate flavours. 
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